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Abstract 

Background: Negative energy balance (NEB), an altered metabolic state, occurs in early postpartum dairy cattle 
when energy demands to support lactation exceed energy intake. During NEB the liver undergoes oxidative stress 
and increased breakdown of fatty acids accompanied by changes in gene expression. It is now known that micro 
RNAs (miRNA) can have a role in mediating such alterations in gene expression through repression or degradation 
of target mRNAs. miRNA expression is known to be altered by metabolism and environmental factors and miRNAs 
are implicated in expression modulation of metabolism related genes. 

Results: miRNA expression was profiled in the liver of moderate yielding dairy cattle under severe NEB (SNEB) and 
mild NEB (MNEB) using the Affymetrix Gene Chip miRNA_2.0 array with 679 probe sets for Bos-taurus miRNAs. Ten 
miRNAs were found to be differentially expressed using the 'samr' statistical package (delta = 0.6) at a q-value FDR 
of < 12%. Five miRNAs including miR-17-5p, miR-31, miR-140, miR-1281 and miR-2885 were validated using RT-qPCR, 
to be up-regulated under SNEB. Liver diseases associated with these miRNAs include non-alcoholic fatty liver (NAFLD) 
and hepatocellular carcinoma (HCC). miR-140 and miR-17-5p are known to show differential expression under oxidative 
stress. A total of 32 down-regulated putative target genes were also identified among 418 differentially expressed 
hepatic genes previously reported for the same animal model. Among these, GPR37 (G protein-coupled receptor 
37), HEYL (hairy/enhancer-of-split related with YRPW motif-like), DNJA1, CD14 (Cluster of differentiation 14) and 
GNS (glucosamine (N-acetyl)-6-sulfatase) are known to be associated with hepatic metabolic disorders. In 
addition miR-140 and miR-2885 have binding sites on the most down-regulated of these genes, FADS2 (Fatty 
acid desaturase 2) which encodes an enzyme critical in lipid biosynthesis. Furthermore, /-//VF3-gamma (Hepatocyte 
nuclear factor 3-gamma), a hepatic transcription factor (TF) that is involved in IGF-1 expression regulation and 
maintenance of glucose homeostasis is a putative target of miR-31. 

Conclusions: This study shows that SNEB affects liver miRNA expression and these miRNAs have putative 
targets in hepatic genes down-regulated under this condition. This study highlights the potential role of miRNAs 
in transcription regulation of hepatic gene expression during SNEB in dairy cattle. 



Background 

Over the past few decades improvements in milk pro- 
duction through genetic selection have been associated 
with a reduction in cow fertility and this decrease in 
fertility has become a major concern for dairy producers 
[1,2]. Reproduction is an energetically expensive process 
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and an altered metabolic state called NEB has been 
established as one of the major physiological causes of 
decreased fertility in high yielding dairy cattle [3-5]. NEB 
is the result of increased energy demands to support lac- 
tation, coupled with lowered feed intake [4,5]. Immune 
related and hepatic functions are known to be effected 
by NEB [6-9] and there is an increased metabolic load 
on the liver to overcome the energy deficit under NEB. 
There is a dramatic increase in hepatic oxidation of fatty 
acids for energy production [9,10] and in addition there 
are extra demands on the liver to increase glycogenesis 
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to meet the glucose requirements of milk production. In 
order to understand the complex metabolic adjustments 
in postpartum dairy cattle liver during NEB, a dairy cattle 
model based on different milking regimes was developed. 
Earlier in a microarray gene expression study 418 hepatic 
genes were reported to be differentially expressed as a re- 
sult of SNEB [9]. These differentially expressed genes have 
roles in lipid metabolism and glycogenic processes, im- 
mune response and the somatotropic axis involved in milk 
production. The regulation in gene expression under 
SNEB, however, is yet to be fully understood. 

It is now known that a class of small RNAs called 
miRNA, about 19 to 25 nucleotides in length, can regulate 
such alterations at the gene expression level. miRNAs 
were first discovered from the round worm C. Elegans al- 
most 3 decades ago [11]. Since their discovery, hundreds 
of miRNAs have been identified across the plant and 
animal kingdoms. It is estimated that whereas only 1-5% 
of genomic transcripts in mammals code miRNAs, up to 
60% of genes are regulated by miRNAs [12-14]. 

The biogenesis of miRNA is a multistep process that 
begins with miRNA gene transcription resulting in pri- 
mary miRNA (pri-miRNA) in the nucleus, followed by 
the generation of around 70nt long stem-loop precursor 
miRNAs (pre-miRNAs) from pri-miRNA. Pre-miRNAs 
are then translocated to the cytoplasm where they are 
trimmed to remove the terminal loop and release -22 nt 
long duplex mature miRNAs containing a guiding strand 
and a passenger strand. The guiding strand assembles 
into cytoplasmic RNA-induced silencing complex (RISC) 
with argonaute protein that guides the complex to their 
complementary mRNA targets [15,16]. The miRNA target- 
mRNA complementary base-pair interaction generally 
occurs between the target site on the 3'UTR of the 
mRNA and the 2 nd to 8 th nucleotides on the 5'UTR of 
the miRNA called the seed region'. Generally, miRNAs 
repress translation through deadenylation of mRNA 
leading to their subsequent degradation or translation 
repression. miRNAs are reported to regulate a wide 
range of biological processes including cell cycle regulation, 
proliferation and differentiation as well as development, 
immune response, carcinogenesis and various metabolic 
processes [17-19]. The effects of metabolites and environ- 
mental factors including hormones, cytokines and nutri- 
ents on miRNAs is well established [20-24]. Alterations in 
cattle liver miRNA expression in response to external ana- 
bolic steroids was reported and miRNAs were suggested as 
potential biomarker for drug abuse in cattle [24]. In 
another study testosterone treatment was reported to alter 
miR-22, miR-690, miR-122, let-7a, miR-30 and let-7d ex- 
pression in female rat liver [25]. Increases in expression of 
miR-155 and miR-132 due to oxidative stress was reported 
in an ALD (alchololic liver disorder) mouse model [26]. 
Oxidative stress of hepatocytes was also reported to alter 



miR-199a-5p expression [27]. Nutritional modulation of 
miRNA expression has been reported in various dietary 
intervention studies of metabolic disorders like obesity, 
diabetes and fatty liver [28-30]. A high fat diet was reported 
to alter the adipose tissue miRNA expressional profile of 
miR-19a, -92a, -92b, -101, -103, -106, -142-5p, and 296 in 
cattle [31]. miRNAs can regulate metabolism and homeo- 
stasis of high energy metabolites as well as insulin signal- 
ling and glucose homeostasis [32-34]. miR-33a and 
miR-33b located within the sterol regulatory element- 
binding proteins (SREBP), key transcription regulators 
of genes involved in cholesterol biosynthesis and up- 
take, regulate cholesterol homeostasis jointly with their 
host genes [35] and have roles in the regulation fatty 
acid metabolism and insulin signaling [36]. In addition, 
miR-122, a liver specific miRNA regulates hepatic fatty 
acid oxidation and fatty acid and cholesterol synthesis 
rate [32,37]. miRNAs are associated with the patho- 
physiology of hepatic metabolic disorders like NAFLD 
and NASH (non-alcoholic steatohepatitis) as well as 
HCC in mouse and human studies [38-41]. The poten- 
tial role miRNAs have to play in NEB in the dairy cow 
remains to be elucidated. 

This study set out to determine a) if hepatic miRNA 
expression was altered as a result of SNEB in a model of 
moderate yielding Holstein Friesian dairy cows in the 
early postpartum period, and b) to integrate hepatic miRNA 
and mRNA expression profiles through prediction of targets 
of these miRNAs among a set of previously reported differ- 
ential expressed hepatic genes under SNEB from the same 
animal model. Elucidation of the expression patterns of 
miRNAs and computational identification of their putative 
target among genes regulated under SNEB will contribute to 
the understanding of the roles of miRNAs in regulating gene 
expression during SNEB in dairy cows. 

Methods 

Animal model 

The animal model employed in this study has been de- 
scribed previously [5,9]. All procedures were carried out 
under license in accordance with the European Commu- 
nity Directive, 86-609-EC. In brief, multiparous Holstein- 
Friesian cows were blocked according to parity, body 
condition score, and previous milk yield, two weeks 
prior to expected calving and within block were ran- 
domly allocated to two treatments; mild NEB (MNEB, 
n = 12) or severe NEB (SNEB, n = 12) groups based on 
different feeding and milking regimes. On day 2 after 
calving, MNEB cows were fed ad libitum grass silage 
with 8 kg per day of a 21% crude protein dairy concen- 
tration and milked once daily; SNEB cows were fed 25 
kg of silage per day silage with 4 kg of crude protein 
per day and milked three times daily. Three times a 
day versus once a day milking was used to advance a 
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state of SNEB in one group by increasing energy with- 
drawal with a concomitant lower need for differences 
in energy intake between groups. The chemical com- 
position of silage and concentrate offered as previously 
described [5]. For sample collection, cows were se- 
lected from each group based on extremes of energy 
balance (MNEB, n = 6; SNEB, n = 6). Cows were slaugh- 
tered approximately 14 days postpartum (MNEB; 13.6 ± 
0.75, range 11-15; SNEB 14.3 ± 0.56, range 13-16) [9]. 
Liver tissues were retrieved within 30 min of slaughter 
and snap frozen at -80°C. One MNEB cow was removed 
retrospectively as she was deemed to be ill during the 
experimental period. 

RNA extraction and quality analysis 

For the purpose of RNA extraction, 100 mg frozen liver 
tissue was directly immersed in 1 ml Trizol (Invitrogen, 
Stockholm, Sweden). A Precellys 24 homogeniser (Bertin 
Technologies, Montigny-le-Bretonneux, France) was used 
to homogenise the tissue (6000 rpm for 20 sec). Hom- 
ogenate was left for 5 min at RT to permit the complete 
dissociation of nucleoprotein complexes followed by cen- 
trifugation at 12,000 g for 10 min at 4°C. Supernatant 
was transferred to a 2 ml tube. 2ul of Pellet paint® NF 
Co-Precipitant (EMD Millipore, Darmstadt, Germany) 
and 200 ul chloroform was added and incubated at RT 
for 5 min. After centrifugation for 15 min at 4°C and 
14000 g the aqueous phase was transferred into a new 
tube and 500 ul isopropanol was added and incubated 
for 5 min on ice. Following centrifugation for 10 min 
at 4°C and 14000 g the supernatant was removed and 
the pellet washed with 500 ul 75% ETOH followed by 
centrifugation for 10 min at 4°C and 14000 g. A 100 ul 
aliquot of nuclease-free water was added to the pellet 
and incubated at RT for 10 min. After vortexing for 
1 min 10 ul of 3 M sodium acetate was added. Next, 
250 ul ice cold ETOH was added and incubated for 
1 hr at -20°C, followed by centrifugation at 4°C for 
10 min and 14000 g. The supernatant was removed 
and washed with 500 ul of 75% EtOH by vortexing 
followed by centrifugation at 4°C for 10 min and 
14000 g. The supernatant was removed and the pellet 
was air-dried for 5 min at RT following which 30ul of 
RNase-free water was added and left for 5 min. RNA 
quality and quantity was assessed using a NanoDrop 
spectrophotometer (ND-1000; Wilmington, DE, USA) 
and Bioanalyzer 2100 (Agilent Technologies Ireland, 
Dublin, Ireland). Mean RIN values were > 8. Total RNA 
was stored at -80°C. 

Gene Chip miRNAs_2.0 array hybridization 

Total RNA from liver samples of 11 cows were hybridized 
to the Affymetrix Gene Chip miRNA_2.0 (Affymetrix, 
Santa Clara, CA, USA) arrays for miRNA expression 



profiling. This array has miRBase vl5 coverage (www. 
mirbase.org) with 15,644 probe sets of 131 organisms 
including 679 probe sets for Bos-taurus miRNAs. Array 
hybridization was carried out at ATLAS Biolabs, Berlin, 
Germany. Briefly, 200 ng total RNA labelled with Flash 
Tag Biotin was hybridized to the miRNA Array in a 
Hybridization Oven 640 (Affymetrix) at 48°C for 18 h. 
The arrays were stained with the Fluidics Station 450 
using fluidics script FS450_0003 (Affymetrix), and then 
scanned on a GeneChip Scanner 3000 7G (Affymetrix, 
Santa Clara, CA, USA). 

Statistical analysis 

All statistical analyses were performed in the R statistical 
computing environment (Version 2.14; http://www.r-project. 
org) with the samr package from the Bioconductor project 
(http://www.bioconductor.org). Data quality was assessed 
with the ArrayQualitymetrix package from Bioconductor 
[42,43]. One eel file from each group was discarded due to 
issues with hybridisation quality. The expresso method of 
the Affy package [44] was used to pre-processes the 
data with quantile normalization and median polish 
summarization. While the chip contains probes com- 
plementary to miRNA from other species only the 
expression intensities of 679 cattle miRNAs were used 
for further analysis. The Significance Analysis of Micro- 
array (SAM) two-class unpaired method implemented in 
the samr Bioconductor package was used to identify 
differentially expressed miRNAs [45]. The samr par- 
ameter delta 0.6 was selected that classified miRNAs 
as differentially expressed at a fold-change > 1.25 and a 
q-value < 12%. 

RT-qPCR validation of differentially expressed miRNAs 

Seven candidate miRNAs differentially expressed on the 
microarray; including miR-17-5p, miR-31, miR-140, miR- 
1281, miR-2885, miR-296, and miR-671 were selected for 
RT-qPCR validation. miRNA expression was carried out 
with TaqMan miRNA RT-qPCR assays according to the 
manufacturers instructions (Applied Biosystems, Dublin, 
Ireland). miRNA-specific reverse transcription was per- 
formed on 10 ng of purified total RNA using the TaqMan 
MicroRNA Reverse Transcription kit according to manu- 
facturers instructions. RT-qPCR reactions were performed 
using 1 ul of cDNA (10 ng/ul) in 9 ul of Taqman universal 
master mix containing TaqMan PCR primers and probes 
on a BioRad CFX96 real time PCR system (Bio-Rad, 
Hemel Hempstead, UK) using the following cycling pa- 
rameters; 95°C for 10 min followed by 40 cycles at 95°C 
for 15 s and 60°C for 1 min. Four biological replicates 
from each of the MNEB and SNEB treatment groups were 
used for RT-qPCR validation. Three reference miRNAs in- 
cluding miR-122, let-7b and RNU6B were tested with the 
software program geNorm version 3.5 [46] for calculating 
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the gene expression stability measure (M value). RNU6B 
was found to be the most stable internal reference miRNA 
across treatments with an M value of 0.7. It was more 
stable on its own than when used in combination with the 
two other reference miRNA. The software package BioRad 
CFX manager was used for correction of the Ct values 
and normalization to RNU6B using the 2" AACt method 
[47]. Corrected Ct values were used to calculate differen- 
tial expression using the PROC t-test (SAS) [48]. 

Prediction of differentially expressed miRNAs targets 
among differentially expressed hepatic genes under SNEB 

A target prediction algorithm for custom data, from the 
TargetScan database website [13] implemented in Perl 
(targetscan_61_context_scores.pl), was used to identify 
target sites of differentially expressed miRNAs in the 
3'UTRs of 418 differentially expressed genes reported 
previously for the same liver tissue [9]. Using this algo- 
rithm, conserved 8mer and 7mer sites matching the 
seed region of each miRNA in the 3'UTRs of differen- 
tially expressed hepatic genes were identified [49]. Non- 
conserved sites and sites with mismatches in the seed 
region that were compensated by conserved 3 ' pairing 
were also included [13]. Predictions were ranked based 
on the context + scores of the sites that represent the 
predicted efficacy of targeting [50]. GO (Gene Ontol- 
ogy) biological processes associations of 3 conserved 
miRNAs (miR-31, miR 17-5p and miR-140) were re- 
trieved using FAME (functional assignment of miRNAs 
via enrichment) software [51]. No data was available 
for miR-1281 and miR-2885 in the FAME database. 
The FAME algorithm makes direct inference of a 
specific miRNA function using enriched subsets of the 
target genes sharing a common biological process or 
pathway. The GO biological processes associated with 
down-regulated putative targets were retrieved from 
the UniProt Gene Ontology Annotation (GOA) data- 
base (http://www.ebi.ac.uk/GOA). 

Results 

Ten out of 679 mature bovine miRNAs represented on 
the miRNA Gene Chip miRNA_2.0 array were found to 
be up-regulated in SNEB cows. Table 1 lists the miRNAs 
that showed differential expression between MNEB and 
SNEB groups. RT-qPCR results confirmed the differen- 
tial expression of 5 out of the 7 miRNAs tested. The 
expression of miR17-5p, -1281, -140, 2885 and -31 were 
consistent between microarray and RT-qPCR whereas 
miR-296 and miR-671 were not significant. Table 2 lists 
the differentially expressed miRNAs validated by RT- 
qPCR. Putative targets of these five miRNAs were identi- 
fied in the set of 418 (230 down- and 188 up-regulated) 
hepatic genes reported previously in the same animal 
model. Among the 418, 67 (35 up-regulated and 32 



Table 1 Differentially expressed post-partum dairy cattle 
miRNAs under SNEB 



miRNA 


miRBase accession 


q-value (%) 


Fold change 


bta-miR-31 


MI0004762 


0.00 


1.79 


bta-miR-1281 


MI0010466 


11.97 


1.71 


bta-miR-2483 


MI001 1545 


11.97 


1.56 


bta-miR-2885 


MI0013058 


0.00 


1.55 


bta-miR-296 


MI0009786 


11.97 


1.52 


bta-miR-2316 


MI001 1336 


11.97 


1.40 


bta-miR-140 


MI0005010 


11.97 


1.33 


bta-miR-17-5p 


MIMAT0003815 


11.97 


1.30 


bta-miR-671 


MI0009887 


11.97 


1.29 


bta-miR-2455 


MI001 1512 


11.97 


1.28 



Differential expression analysis was done for microarray expression intensities 
of 679 bovine miRNAs with SAM (significance analysis of microarrays) method 
at delta parameter of 0.6 and q FDR < 12%. 



down-regulated) genes were found to have target sites in 
the 3'UTRs of the up-regulated miRNAs (Table 3). 

All of the five differentially expressed miRNAs vali- 
dated by RT-qPCR in this study were up-regulated how- 
ever there was no significant enrichment among the 
down-regulated putative targets genes using the Fisher 
test for enrichment analysis. Overall miR-140 was found 
to have the most putative targets among differentially 
expressed genes under SNEB while miR-17-5p had the 
second highest number of targets. As the consensus is 
that up-regulated miRNAs result for the most part in 
down-regulation of their target genes, we further se- 
lected only down-regulated genes as putative affected 
targets for subsequent analysis to investigate the role 
of differentially expressed miRNAs under SNEB. As 
miRNAs are thought not to switch off their targets 
completely but rather fine-tune their expression, all 
down-regulated genes were included irrespective of 
fold-change. Figure 1 presents the subset of the puta- 
tive down-regulated targets for these up-regulated 
miRNAs. miR-140 had the highest number of putative 
down-regulated target genes followed by miR-17-5p. 



Table 2 Differentially expressed postpartum dairy cattle 
hepatic miRNAs under SNEB 



miRNA 


Fold change 


p value 


bta-miR-31 


4.11 


0.033 


bta-miR-1281 


3.10 


0.042 


bta-miR-2885 


2.91 


0.006 


bta-miR-17-5p 


4.62 


0.008 


bta-miR-140 


3.90 


0.030 


bta-miR-296 


0.85 


0.490 


bta-miR-671 


0.92 


0.930 



Differential expression analysis was carried out with PROC t-test (SAS) after 
normalization of Ct-values to reference gene RNU6B. 
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Table 3 Putative target genes of up-regulated hepatic miRNAs differentially expressed under SNEB 



Down-regulated putative targets Up-regulated putative targets 



Symbol 


Gene name 


Symbol 


Gene name 


ALDH1A1 


aldehyde dehydrogenase 1 family, member A1 


AADAT 


aminoadipate aminotransferase 


BTG1 


B-cell translocation gene 1, anti-proliferative 


ACP2 


Acid phosphatase 2, lysosomal 


CCL19 


chemokine (C-C motif) ligand 19 


ADTRP 


Androgen-Dependent TFPI-Regulating Protein 


GDI 4 


CD14 molecule 


ARG1 


Arginase, liver 


CIB1 


Calcium and integrin binding 1 (calmyrin) 


BIN1 


Bridging integrator 1 


DNAJA1 


DnaJ (Hsp40) homolog, subfamily A, member 1 


BNIP3L 


BCL2/adenovirus E1B 19 kDa interacting protein 3-like 


DSG1 


desmoglein 1 


CMPK1 


cytidine monophosphate (UMP-CMP) kinase 1, cytosolic 


EPB41L5 


Erythrocyte membrane protein band 4.1 like 5 


CNN1 


Calponin 1, Basic, Smooth Muscle 


ERRFI1 


FRRR rarpntnr fppHh^rk inhihitnr 1 

l_r\L)L) ICL.C|-JIUI ICCUUC1LN II II IIUIIUI I 




Ldl UUAy |Jc|JUUcl be 


FADS 2 


Fsttx/ ^riH Hp^^ti ir^^p 7 
i any u^iu ucjului ojc z. 


CPT1B 


r^rnitinp n^lmitnvltr^n^fpr^^p 1R (mi kHp^ 
iwCiiiiiuiic pa 1 1 1 1 1 L^y i li a i i ji ci ajc i u \jiiujv_ic/ 


F0XA3 


Forkhead box A3 


DENND2D 


DENN/MADD domain containing 2D 


FXR1 


Fragile X mental retardation, autosomal homolog 1 


DPYD 


dihydropyrimidine dehydrogenase 


GKAP1 


G kinase anchoring protein 1 


HDC 


histidine decarboxylase 


GNS 


glucosamine (N-acetyl)-6-sulfatase 


HMGN4 


High mobility group nucleosomal binding domain 4 


GPBP1 


GC-rich promoter binding protein 1 


HSF2BP 


Heat shock transcription factor 2 binding protein 


GPR37 


G protein-coupled receptor 37 

(pnHothplin rprpntnr fvnp R-likp^ 


ILIA 


Interleukin 1, alpha 


HEYL 


hairy/enhancer-of-split related with YRPW motif-like 


LPCAT3 


I v^nnhn^nhatiHvlrhnlinp arvllrarKfpr^KP 3 


IP6K2 


inositol hexakisphosphate kinase 2 


MBL2 


Mannose-binding lectin (protein C) 2, soluble 
(opsonic defect) 


KIAA1 191 


KIAA1 191 


MINA 


MYC induced nuclear antigen 


MIEN1 


Migration and invasion enhancer 1 Bottom of Form 


NHEJ1 


Nonhomologous end-joining factor 1 


MIS 12 


MIND kinetochore complex component, homolog (S. pombe) 


OAT 


ornithine aminotransferase (gyrate atrophy) 


NR2F1 


nuclear receptor subfamily 2, group F, member 1 


PDRG1 


p53 and DNA-damage regulated 1 


PRSS35 


protease, serine, 35 


P0LR3GL 


polymerase (RNA) III polypeptide G (32kD)-like 


PTPRR 


protein tyrosine phosphatase, receptor type, R 


RCAN3 


RCAN family member 3 


RNF 128 


ring finger protein 128 


REEP5 


Receptor accessory protein 5 



Putative binding sites were identified in the 3'UTR of differentially expressed hepatic genes under SNEB with the Targetscan algorithm. 



miR-31 was the most up-regulated miRNA under SNEB 
while miR-1281 had all but one of its putative target genes 
down-regulated under SNEB. miR-2885 has FADS2 as its 
only putative down- regulated target under SNEB and 
FADS2 is also a putative target of miR-140. DSG1 (desmo- 
glein 1) and FXR1 (fragile X mental retardation, autosomal 
homolog 1) are both shared targets of miR-17-5p and 
miR-31 while GNS (glucosamine-N-acetyl-6-sulfatase) is a 
shared target of mir-140 and miR-31 (Figure 2). Table 4 
lists the bio-types of putative targets down-regulated 
under SNEB. The GO biological processes associated with 
down-regulated putative targets are given in Additional 
file 1: Table SI. The overall GO biological functions asso- 
ciated with three conserved miRNAs that are relevant in 
terms of SNEB are given in Additional file 1: Table S2. 
miR-31 is associated with the GO biological processes of 
cell growth, amino acid transport and response to stress. 



miR-17-5p is associated with GO biological processes 
include regulation of cell cycle and protein metabolism. 
miR-140 is associated with GO biological processes of cell 
proliferation, polysaccharide metabolism and signal trans- 
duction, while both miR-31 and miR-140 are associated 
with response to nutritional level. 

Discussion 

Dairy cows enter a state of NEB postpartum when the 
energy supply is prioritised towards the mammary 
tissues for milk production while the energy intake is 
reduced. The principal physiological response to the 
energy deficit is mobilization of non-esterified fatty acids 
(NEFA) from adipose tissue. NEFAs are broken down in 
the liver to release energy and any incomplete break- 
down results in production of BHB (p-hydroxybutyrate) 
or TAGs (triacylglycerides). Accumulation of TAGs, BHB 
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miR-140 


FC 


mirl7-5p | FC 


bta-miR-31 


FC 


mir2885 


FC 


miR-1281 


FC 


FADS 2 


-19.60 


DSG1 [ -2.10 


DSG1 


-2.10 


FADS 2 


-19.60 


HEYL 


-2.27 


CCL19 


-1.50 


GPR37 


-1.84 


FOXA3 


-1.48 


RNF4 


1.19 


SLC1A4 


-2.20 


BTG1 


-1.44 


SYBU 


-1.72 


PTPRR 


LPCAT3 


1.25 


IP6K2 


-1.43 


TSC22D1 


-1.41 


PRSS35 


-1.50 


MIS 12 


-1.42 


BIN1 


1.47 


NR2F1 


-1.19 


KIAA1191 


-1.32 


TTR 


-1.32 


CD14 


-1.39 


TUFT1 


2.05 


POLR3GL 


1.15 


DNAJA1 


-1.31 


FXR1 


-1.26 


FXR1 


-1.26 


SCARB 1 


2.62 






MIEN1 


-1.28 


RNF128 


-1.24 


NHEJ1 


1.27 


SLC25A34 


3.29 






EPB41L5 


-1.25 


TPMT 


-1.24 


ARG1 


1.32 










ERRFI1 


-1.25 


GKAP1 


-1.21 


SEC61A2 


1.41 










CIB1 


-1.15 


ALDH1A1 


-1.12 


OAT 


1.59 










GNS 


-1.12 


GNS 


-1.12 


HSF2BP 


1.80| 








STAMBP 


-1.09 


CMPK1 


1.15 


SLC25A48 


1.93 










TMX1 


-1.09 


HMGN4 


1.17 














GPBP1 


-1.08 


RNF4 


1.20 














DPYD 


1.12 


BNIP3L 


1.20 














MINA 


1.13 


SLC27A4 


1.25 














SNAPIN 


1.15 


CNN1 
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Figure 1 Predicted miRNA gene targets with fold-change under SNEB. 



and NEFAs results in the onset of oxidative stress in the 
liver that can lead to fatty liver or lipidosis. In addition 
components of the somatotrophic axis like IGF (insulin 
like growth factor) and GH (growth hormones) that have 
roles in the control of milk production are also affected 
by SNEB [8-10]. Reduced IGF1 expression is also nega- 
tively correlated with fertility in dairy cows [52]. Previ- 
ous nutritional supplementation of animals in NEB with 
lipogenic and glycogenic nutrients partly succeeded in 
compensating for energy partitioning between mammary 



tissues and the rest of the body [4,53]. However, the 
regulation of such alterations is yet to be explored. In 
this study, five miRNAs were altered in response to 
SNEB. Two of these miRNAs have been reported to be 
altered under oxidative stress in humans [41,54] and 
miR-17-5p has been associated with oxidative stress in 
HCC [41]. In addition alterations in miR-140 expression 
under chemical and oxidative stress has also been re- 
ported [54]. Three of the miRNAs, miR- 17, miR-31 and 
miR-140, found to be up-regulated under SNEB have 
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Table 4 Intracellular location and bio-types of down-regulated putative gene targets of up-regulated miRNAs 



Symbol 


Entrez gene name 


Location 


Type 


|/| A A 1 1 Q1 
i\Iaa/A I I y I 


klA A1 1 Q1 
MAM I I y I 


Cytoplasm 


vjii lei 


K/IIFM 1 
IVIILlN I 


IVIiy IdllUI 1 dllU INVdblON clllldllLcf I 


L.y LCjpidSI 1 1 


Ull lcl 


IV /IK 1 9 
IVID I Z 


MIS 12, MIND kin6tochor6 complex component, homolog (S. pomb6) 


Nucleus 


Other 

uiner 


MR9F1 
INnZr I 


Nuclear receptor subfamily 2, group F, m6mb6r 1 


Nucleus 


Ligand-dependent nuclear receptor 


DDCCoc 
r njjjj 


rlULcdbc, bcllllc, jj 


LXLI dLcllUldl bpdLc 


pcpilUdbc 


PTPRR 


rlULclll LyiUblllc pi lUbpi IdLdbc, IcLcpLUI Lyptr, n 


rldbllld 1 1 lcl I lUldl lc 


pi lUbpi Id Ldbc 


RMF 1 9R 
nl\l r I zo 


nlliy IIIILjcl piULclll I Zo, Co UUiyUlllll piUlclll IILjdbc 


^_y LUpidbl 1 1 


li izyi i ic 


SLC1 A4 


Qpih ito warrior famik/ 1 fpili itamato/noi itral aminn ar~iH trancpiortor^ momKor A. 
OUIULc Ldlllcl Idlllliy I vy lULdl I IdLcV I IcULI dl dllllllU dLIU LI dl lo|JUI Lcl ), lllclllUcI H 


r Idol 1 Id 1 1 lcl 1 IUI dl Ic 


1 1 dl lo|JUI Lcl 


_> I AlVIDr 


0 1 AIVI UlllUlliy piULclll 


Mi i r~\ i c 
INULIcUS 


F r~\ ~7\ / rv~\ o 
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1 1 rt Lr\r^\A/r\ 
U 1 IKI IUVVI 1 


Plthpr 
ULI lcl 


TMX1 


Thioredoxin-related transmembrane protein 1 


Cytoplasm 


Enzyme 


TPMT 


Thiopurine S-methyltransferase 


Cytoplasm 


Enzyme 




Transthyretin! 


Cytoplasm 


Transporter 


ALDH1A1 


Aldehyde dehydrogenase 1 family, member A1 


Cytoplasm 


Enzyme 


BTG1 


B-cell Translocation gene 1, anti-proliferative 


Nucleus 


Transcription Regulator 


CCL19 


chemokine (C-C motif) ligand 19 


Extracellular space 


Cytokine 


CD14 


CD14 molecule 


Plasma Membrane 


Transmembrane receptor 


CIB1 


Calcium and integrin binding 1 (calmyrin) 


Nucleus 


Other 


DNAJA1 


DnaJ (Hsp40) homolog, subfamily A, member 1 


Nucleus 


Other 


DSG1 


Desmoglein 1 


Plasma membrane 


Membrane protein 


EPB4 1 L5 


Erythrocyte membrane protein band 4.1 like 5 


Plasma membrane 


Membrane protein 


ERRFI1 


ERBB receptor feedback inhibitor 1 


Cytoplasm 


Other 


FADS2 


Fatty acid desaturase 2 


Plasma membrane 


Enzyme 


F0XA3 


Forkhead box A3 


Nucleus 


Transcription regulator 


FXR1 


Fragile X mental retardation, autosomal homolog 1 


Cytoplasm 


Other 


GKAP1 


G kinase anchoring protein 1 


Cytoplasm 


Other 


GNS 


Glucosamine (N-acetyl)-6-sulfatase 


Cytoplasm 


Enzyme 


GPBP1 


GC-rich promoter binding protein 1 


Nucleus 


Transcription regulator 


GPR37 


G protein-coupled receptor 37 (endothelin receptor type B-like) 


Plasma membrane 


G-protein coupled receptor 


HEYL 


Hairy/enhancer-of-split related with YRPW motif-like 


Nucleus 


Transcription regulator 


IP6K2 


Inositol hexakisphosphate kinase 2 


Cytoplasm 


Kinase 
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been previously associated with hepatic disorders [55-58] 
and miR-140 has been associated with NAFLD [40,57]. A 
regulatory role for miR-17-5p has been reported in HCC 
[55,56] while miR-31 has also been associated with both 
NAFLD and HCC [58-60]. 

A number of the putative down- regulated targets of 
the differentially expressed miRNAs in this study are 
known to have roles in hepatic disorders. These include 
the miR-17-5p target GPR37 (G protein-coupled recep- 
tor 37) a member of the G-protein coupled receptor- 1 
family associated with NAFLD [61]. The miR-1281 tar- 
get HEYL (hairy/enhancer-of-split related with YRPW 
motif-like) is associated with HCC [62,63]. The miR-140 
target DNJA1 that has a role in lipid intake is associated 
with the oxidative stress response of high fat induced 
NAFLD [28,30]. The miR-31 target, CD 14 (Cluster of 
differentiation 14) involved in transmission and release 
of polysaccharides, is also implicated in ALD, hepatic 
cholestasis and hepatic fibrosis [64,65]. The miR-17-5p 
target GNS (glucosamine (N-acetyl)-6-sulfatase) that 
metabolises glucosamine sulphate has been implicated in 
hepatocellular dysplasia, cirrhosis of the liver and hep- 
atic fibrosis [66]. Moreover, the miR-17-5p target TTR, a 
carrier protein for transport of lipid soluble vitamins 
from the liver to the circulation is also reduced under 
SNEB. Reduced TTR expression has been reported in a 
study of periparturient period Holstein and Jersey dairy 
cows in NEB [67]. In addition we were able to postulate 
a role for miR-31 in SNEB through its putative target 
hepatic transcription factor FOXA3 also known as 
HNF3-gamma. Bovine IGF-I has binding sites for HNF3- 
gamma in its promoter region [68]. While another re- 
cent human study reported target sites of miR-31 on the 
3'UTR oilGF-1 [69]. Mf-gamma is also known to be 
involved in glucose homeostasis in hepatocytes through 
regulation of GLUT2 (glucose transporter 2) [70,71]. 
However, the most interesting putative target of all was 
the gene encoding the lipogenesis enzyme FADS2 which 
is critical for long-chain polyunsaturated fatty acids 
biosynthesis [72]. In this study we found binding sites 
for two up-regulated miRNAs miR-140 and miR-2885 
in the 3'UTR of FADS2, whereas miR-140 has been 
previously implicated in the regulation of FADS2 
homolog FADS1 [73]. 

Conclusion 

This study suggests a role for hepatic miRNAs in lipid 
and energy metabolism through the identification of a 
subset of their putative targets associated with important 
metabolic processes. Moreover hepatic miRNAs associ- 
ated with genes involved in the somatotropic axis can be 
a possible link between reduced reproductive perform- 
ance and SNEB in the early postpartum dairy cow. Fur- 
ther direct functional studies of selected miRNA-mRNA 



putative pairs in the liver of such models will help to im- 
prove our understanding of the molecular mechanisms 
and pathways involved in SNEB in postpartum dairy 
cattle. 
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